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This paper discusses a method to integrate a wind turbine generating system 
(WTGS) into a three-phase unbalanced distribution system load flow 
(DSLF) analysis. The proposed method is based on the single-phase 
multiple-node model. In the present work, the single-phase multiple-node 
model is extended to a three-phase multiple-node model to facilitate the load 
flow analysis of a three-phase unbalanced power system network. The 
multiple-node model (i.e., three-node model) will only modify the load flow 
analysis by introducing two lines and two load buses to the distribution 
system network where the WTGS is installed. Thus, a standard three-phase 
load flow program can be employed to compute the unknown quantities in 
the DSLF problem formulation. The proposed method is verified by 
incorporating the model into the load flow analysis of three-phase 
distribution networks. The investigation uses two representative distribution 
networks (i.e., 19-bus and 25-bus networks). The results of the study 


confirm the validity of the proposed method. 
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1. INTRODUCTION 

In general, the electric power distribution system is highly unbalanced. The system unbalance is 
mainly caused by: i) uneven loads on each phase; and/or ii) untransposed distribution lines. Due to the 
unbalance, the single-phase approach can no longer be used to analyze such a system, and a three-phase 
method should be used instead [1], [2]. Moreover, the integration of distributed generations with renewable 
energy sources in the electric power distribution system has also increased recently. It has been 
acknowledged that primary energy sources for distributed renewable energy generation include: wind power, 
solar photovoltaic, and hydropower. However, wind-driven electric power generation systems are more 
established in technology development among the three energy sources and have broader applications 
[3]-[5]. Since the penetration of wind turbine generating system (WTGS) in electric power distribution 
systems has been increasing significantly, the assessment of the system’s steady-state performance will be 
more difficult. Generally, the system steady-state performance assessment is carried out based on the results 
of load flow analysis. The first and probably most critical step in the analysis is to model all system 
components. Several attractive techniques for integrating WTGS in steady-state load flow analysis have been 
proposed in [6]-[23]. 

An exciting method to incorporate the WTGS model into a standard load flow program has been 
proposed in [6]-[8]. The inclusion of the WTGS adds two buses (nodes), two series elements, one shunt 
element, and one load on the existing power system network. A mathematical model of WTGS has been 
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developed in [9]-[11]. This mathematical equation (together with the existing load flow equation) is then 
solved iteratively. Various steady-state load flow models of fixed-speed WTGS have been proposed in [12], 
[13]. The models are developed based on the equivalent circuits of the induction generator. These 
mathematical models are combined with the steady-state models of the system without WTGS and then 
solved simultaneously. Some methods for integrating variable speed WTGS have been proposed in 
[14]-[20]. The WTGS models in [14]-[20], are also derived based on the steady-state equivalent circuits of 
the generator used by the WTGS. Steady-state three-phase models of WTGS have been developed in 
[21]-[23]. The models in [21]-[23] can be incorporated and applied in the load flow analysis of unbalanced 
electric power systems. 

In the methods proposed in [16]—[20] the power system has been assumed to be balanced, and the 
single-phase load flow has been used in the analysis. Since, as previously mentioned, the electric power 
distribution system is generally unbalanced, and the single-phase load flow approach is no longer adequate to 
analyze such a system. Although the proposed methods in [21]—[23] can be used to integrate the WTGS into 
an unbalanced system’s load flow analysis, a standard three-phase load flow program can not be employed to 
obtain the solution. The program has to be modified to include the WTGS mathematical model in the 
calculation. 

Therefore, against the above background, the objective of the present paper is to develop the WTGS 
model that can facilitate the load flow analysis in a three-phase unbalanced power system network without 
the need to modify the standard load flow program. The present work proposes the extension of the 
single-phase multiple-node model [6]-[8] to a three-phase multiple-node model. The model is incorporated 
into the load flow analysis of a three-phase unbalanced distribution network. A standard three-phase load 
flow program can then be employed to compute the unknown quantities in the distribution system load flow 
(DSLF) problem formulation. The rest of the paper is organized as follows: section 2 addresses the single- 
phase multiple-node model of WTGS. The extension of the single-phase multiple-node model to the three- 
phase multiple-node model is given in section 3. Section 4 discusses study cases where validation of the 
proposed method is presented. Finally, some important conclusions of the paper are pointed out in section 5. 


2. MULTIPLE-NODE MODEL OF WTGS 

Figures 1 and 2 show induction generator equivalent circuits of WTGS connected to an electric 
power system at bus k [6]—[10]. It is to be noted that both of the equivalent circuits in Figures 1 and 2 are 
valid and can be used to represent the WTGS induction generator. In Figure 1, Pm is turbine mechanical 
power input, S, is WTGS electric power output, Vs and Vp are stator and rotor voltages, Js and Jp are stator 
and rotor currents, and s is induction generator slip. Also, in Figure 2, Zs and Zr are stator and rotor circuit 
impedances, and Yy is magnetic core circuit admittance. These impedances and admittance are determined 
based on the induction generator resistances/reactances and have the following forms: 


Zs = Rs +jXs (1) 
Zr =Rpt jXr (2) 
Yu = (Re + jXm)/jJRcXm (3) 


where Rs/Xs is resistance/reactance of stator circuit, Rr/Xpr is resistance/reactance of rotor circuit, and R/X» is 
resistance/reactance of magnetic core circuit. 


k fs Rs+jXs RrtjXr Ik 

11 | 
f 

Se=PetjQe 


R- Xm  Re(l-s)/s 


Figure 1. Equivalent circuit of induction generator (per-phase) 
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Figure 2. Equivalent circuit of induction generator in terms of impedances (per-phase) 


In Figure 1, the turbine mechanical power input is represented by the active power consumed by the 
variable resistance (Rr(J-s)/s). By looking at Figure 1, this mechanical power can be expressed as (4): 


Pa = -Rr — lglg (4) 


Based on (4), an equivalent constant power source can therefore be used to represent the variable 
resistance, as shown in Figure 2. Thus, (4) can be eliminated from the load flow equations, and slip s can also 
be excluded from the load flow calculations. It should be noted that the value of turbine mechanical power is 
usually available or known since the turbine manufacturer generally provides the turbine mechanical power 
curve. Based on the power curve, the turbine mechanical powers, for any wind speed values, can then be 
determined. 

Based on the above discussion, the WTGS induction generator model can be converted into a 
three-node model in the load flow analysis [6]-[8]. Figure 3 shows the three-node model of the WTGS 
connected to bus k of an electric power system. In Figure 3, it has been assumed that the power system has n 
buses and m branches. The three-node model application will modify the load flow analysis by adding two 
lines and two load buses to the original system network where the WTGS is installed. Thus, a standard load 
flow program can be employed to compute the unknown quantities in the DSLF problem formulation. 
Modification to the source code program does not need to be carried out. What needs to be done is to change 
only the program data input due to the additional lines and buses. 


n+l Zr n+2 


POWER 
SYSTEM 


Yu Pp 


Figure 3. Multiple-node model of WTGS 


Detail of the branch and bus data added to the original system network are shown in Tables 1 and 2. 
The parameters in Tables 1 and 2 are the additional data input to the load flow program due to the WTGS 
installation at bus k. It can be seen from Table 1 that the turbine mechanical power is treated as a load with 
negative power (-P,,). The value of Pm, as mentioned earlier, can be determined based on the power curve 
provided by the turbine manufacturer. It should also be noted that Zs, Zr, and Yy are calculated using (1). 


Table 1. Addition of bus data Table 2. Addition of branch data 
Bus Power Power Shunt Line Bus p toq Series impedance Half line shunt admitance 
generation demand admittance (Z) (Y2) 
n+l 0 0 Yu m+1 kton+l1 Zs 0 
n+2 0 -Pm 0 m+2 n+l ton+2 Zr 0 
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3. MULTIPLE-NODE MODEL OF WTGS FOR THREE-PHASE DSLF ANALYSIS 

Figure 4 shows the application of the three-node model in a three-phase electrical power system 
network. Figure 4 is obtained based on Figure 3, where the impedances of the induction generator (i.e., Zs, 
Zr, and Yy) have been represented by 3x3 matrices to facilitate the analysis in a three-phase network. Detail 
of the branch and bus parameters added to the power system due to WTGS installation are presented in 
Tables 3 and 4. In Tables 3 and 4, the elements of 3x3 series impedance and shunt admittance matrices are 
calculated as (5)-(7): 


ZE = RP + jX (5) 
ZH? = REP + jXH 6) 
YHP = (REY + XEROX ) 


where u=a, b, c and v=a, b, c. Also, the elements of turbine mechanical power vector can be determined 
using: 


P4 = P pe se (8) 


POWER 
SYSTEM 


Figure 4. Multiple-node model for three-phase power flow 


Table 3. Addition of bus data for three-phase system 
Bus _ Power generation 
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0 
0 


Shunt admittance 
f ya | 


Power demand 
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Table 4. Addition of branch data for three-phase system 
Line Busptog Series impedance (Z) _ Half line shunt admittance (Ysh/2) 


4. 


m+1 kton+1 2g? 


m+2 n+lton+2 


zg 
zg: 
Ze 
ZR 
ZR 
Zé 


ooo ooo 


0 
o| 
0 


0 
o| 
0 


CASE STUDY 


4.1. Case study 1 

The first case study is based on 19-bus network. A single-line diagram of the network is shown in 
Figure 5. It is an unbalanced network with a system voltage of 11 kV [23]. The system has a total three-phase 
load of 515.94 kW and 177.27 KVAR (phase-a: 176.33 kW and 61.23 kKVAR; phase-b: 166.34 kW and 
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56.34 kVAR; phase-c: 173.27 kW and 59.70 kVAR). Detail of the system data including the WTGS data can 
be found in [23]. Results of the load flow analysis for the system are shown in Tables 5-7. In Tables 5-7, the 


values in pu are ona | MVA base. It is to be noted that the results in Tables 5-7 are in excellent agreement 


with the results in [23], which confirm the correctness of the proposed method. 


Figure 5. 19-bus network 


Table 5. 19-bus network WTGS voltages (in pu) 


Pn (pu) Phase-a Phase-b Phase-c 
0.03 0.9178721.09° 0.91590.2-118.91° 0.91671.2121.05° 
0.06 0.92543 71.299 0.92352-118.71° 0.92431121.25° 
0.09 0.9328621.49° 0.931002-118.51° 0.93176.2121.45° 
0.12 0.940151.68° 0.938352Z-118.32° 0.93908.2121.64° 
0.15 0.9473221.87° 0.94556.Z-118.13° 0.94627.2121.84° 
0.18 0.9543722.06° 0.95265.2-117.94° 0.95335.2122.03° 
0.21 0.9613122.25° 0.95963.Z-117.75° 0.96030.2122.21° 
0.24 0.9681322.43° 0.966502-117.57° 0.96715.2122.40° 
0.27 0.97485 22.629 0.973262Z-117.38° 0.97389.2122.58° 
0.30 0.9814722.80° _0.97991Z-117.20° _0.98053.2122.76° 


Table 6. 19-bus network WTGS power outputs (in pu) 


Pn (pu) Phase-a Phase-b Phase-c 
0.03 0.00831-j0.03369 0.00832-j0.03355 0.00832-j0.03360 
0.06 0.01828-j0.03425 0.01829-j0.03411 0.01829-j0.03417 
0.09 0.02825-j0.03482 0.02826-j0.03468 0.02826-j0.03474 
0.12 0.03822-j0.03538 0.03823-j0.03525 0.03823-j0.03530 
0.15 0.04819-j0.03594 0.04820-j0.03581 0.04820-j0.03586 
0.18 0.05816-j0.03650 0.05817-j0.03637 0.05816-j0.03643 
0.21 0.06813-j0.03706 0.06814-j0.03694 0.06813-j0.03699 
0.24 0.07810-j0.03762 0.07810-j0.03750 0.07810-j0.03755 
0.27 0.08806-j0.03818 0.08807-j0.03806 0.08807-j0.03810 
0.30 0.09803-j0.03874 _0.09804-j0.03861 0.09803-j0.03866 


Table 7. 19-bus network substation power outputs (in pu) 


Pn (pu) Phase-a Phase-b Phase-c 
0.03 0.17969+j0.09995 0.16982+j0.09493 0.17684+j0.09843 
0.06 0.16840+j0.09995 0.15849+j0.09492 0.16552+j0.09841 
0.09 0.15730+j0.10003 0.14735+j0.09499 0.15440+j0.09848 
0.12 0.14638+j0.10019 0.13640+j0.09514 0.14346+j0.09862 
0.15 0.13564+j0.10042 0.12562+j0.09536 0.13270+j0.09884 
0.18 0.12506+j0.10072 0.11502+j0.09566 0.12210+j0.09912 
0.21 0.11465+j0.10108 0.10458+j0.09602 0.11167+j0.09948 
0.24 0.10439+j0.10152 0.09429+j0.09645 0.10140+j0.09990 
0.27 0.09428+j0.10201 0.08415+j0.09694 0.09127+j0.10039 
0.30 0.08431+j0.10257 0.07416+j0.09749 0.08128+j0.10093 


4.2. Case study 2 


The second case study is based on 25-bus network. A single-line diagram of the network is given in 
Figure 6. It is also an unbalanced network with a system voltage of 4.16 kV [24], [25]. The system has a total 
three-phase load of 11,145 kW and 2,884 KVAR (phase-a: 1,298.3 kW and 957 kVAR;; phase-b: 3.919,9 kW 
and 967 kVAR; phase-c: 1,303.3 kW and 960 KVAR). Al of the system data can be found in [24], [25]. A 
WTGS with a power rating of 3 MW is assumed to be connected to bus 25. Data of the WTGS (i.e., 
induction generator impedances) are shown in Table 8. 
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Figure 6. 25-bus network 


Tables 9-12 show the load flow analysis results for the unbalanced 25-bus distribution system. In 
Tables 9-12, the values in pu are also on a 1 MVA base. Since the system is unbalanced, the electrical 
quantities of the three phases are not identical. The results in Tables 9-12 are also given in the form of 
graphics as shown in Figures 7-9. Figure 7 indicates that the system voltage profile improves as the amount 
of turbine mechanical power is raised. This improvement is due to increased WTGS or wind power plant 
(WPP) output. This WTGS output increment implies that the power injection at the WTGS bus also 
increases, reducing the line losses and improving the system voltage profile. 


Table 8. WTGS data 


Parameter 


Resistance/reactance value (Ohm) 


Stator circuit 


Rotor circuit 


Core circuit 


Self: Rs“=Rs”’=Rs=0.057685; Xs“4=X "=X 5°=0.288427 


Mutual: R”=R =R” =0; X” =X =X =0 


Self: RHR? =RR=0.057685; Xr” =X RP =X" =0.288427 


Mutual: Re’=Re“=Re=0; XR” =X =X" =O 
Self: Ri“ =R =R =576.85; Xi =X,” =K yh =28.84 


Mutual: (RPV =REJ REN = 0; Knit) Kyl =X, 10 


Table 9. 25-bus network WTGS voltages (in pu) 


Pn (pu) Phase-a Phase-b Phase-c 
0.3 0.92143 Z 0.26° 0.92619Z-119.54° 0.934732 119.90° 
0.6 0.92937 Z 0.74? 0.93421Z-119.10° 0.941362 120.35° 
0.9 0.93698 Z 1.22? 0.94185Z-118.65° 0.947672 120.80° 
1:2 0.94402 Z 1.70° 0.949147-118.21° 0.95365% 121.25° 
1.5 0.95078 22.199 0.956102Z-117.76° 0.960332 121.70° 
1.8 0.95720 Z 2.67? 0.96275Z-117.31° 0.96673 Z 122.16° 
2.1 0.96329 Z 3.16° 0.96910Z-116.86° 0.972852 122.61° 
2.4 0.96908 Z 3.65°  0.97517Z-116.41° 0.979702 123.07° 
27. 0.97456 Z 4.14° 0.980972Z-115.96° 0.985302 123.53° 
3.0 0.97976 Z 4.63? _0.98651Z-115.50° _0.99064.2 124.00° 


Table 10. 25-bus network WTGS power outputs (in pu) 


Pn (pu) Phase-a Phase-b Phase-c 
0.3 0.09110-j0.16958 0.09102-j0.17131 0.09086-j0.17444 
0.6 0.19025-j0.17626 0.19017-j0.17800 0.19005-j0.18058 
0.9 0.28897-j0.18495 0.28890-j0.18669 0.28881-j0.18874 
1.2 0.38729-j0.19556 0.38723-j0.19729 0.38717-j0.19882 
1.5 0.48521-j0.20800 0.48517-j0.20970 0.48514-j0.21075 
1.8 0.58276-j0.22219 0.58274-j0.22386 0.58274-j0.22446 
2.1 0.67996-j0.23810 0.67997-j0.23970 —0.67997-j0.23991 
2.4 0.77680-j0.25565 — 0.77684-j0.25717 — 0.77684-j0.25755 
2.7 0.87330-j0.27482 0.87339-j0.27622 0.87336-j0.27685 
3.0 0.96946-j0.29557 _0.96961-j0.29682 _0.96955-j0.29728 
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Table 11. 25-bus network substation power outputs (in pu) 


Pn (pu) 


Phase-a 


Phase-b 


Phase-c 


0.3 


1.28992+j1.22132 
1.18524+j1.22279 
1.08312+j1.22438 
0.98340+j1.23188 
0.88496+j1.24309 
0.79071+j1.25786 
0.69755+j1.27605 
0.60641+j1.29756 
0.51723+j1.32228 
0.42996+j1.35016 


1.31299+j1.22448 
1.20819+j1.22568 
1.10594+j1.22878 
1.00609+j1.23658 
0.90752+j1.24789 
0.81313+j1.26256 
0.71983+j1.28044 
0.62855+j1.30144 
0.53922+j1.32544 
0.45180+j1.35239 


1.27542+j1.22218 
1.17266+j1.22301 
1.07202+j1.22575 
0.97341+j1.23323 
0.87476+j1.24429 
0.78198+j1.25880 
0.68901+j1.27666 
0.59781+j1.29779 
0.50834+j1.32210 
0.42056+j1.34954 


Oo 699 


Table 12. 25-bus network WTGS/substation total outputs and line losses (in pu) 


Pm (pu) WPP Substation line losses 
0.3 0.27298-j0.51533 3.87833+j3.66798 0.23141+j0.26865 
0.6 0.57047-j0.53484 3.56609+j3.67148 0.21666+j0.25264 
0.9 0.86668-j0.56038 3.26108+j3.67891 0.20786+j0.23453 
1.2 1.16169-j0.59167 2.96290+j3.70169 0.20469+j0.22602 
1.5 1.45552-j0.62845 2.66724+j3.73528 0.20286+j0.22283 
1.8 1.74824-j0.67051 2.38582+j3.77922 0.21416+j0.22471 
2.1 2.03990-j0.71771 2.10639+j3.83316 0.22639+j0.23145 
2.4 2.33048-j0.77037 1.83277+j3.89678 0.24335+j0.24241 
2.1 2.62005-j0.82789 — 1.56479+j3.96983 0.26494+j0.25794 
3.0 2.90862-j0.88967 1.30232+j4.05208 0.29104+j0.27841 
1— — — r — 
0.99 
oO 
goe 
S 0.97 
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e 0.95 
a 
$ 0.94 
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Figure 7. Variation of WPP voltages 


Figures 8(a) and (b) show the variations of WTGS power outputs against turbine mechanical power 
inputs. Figure 8(a) show that with the increase in turbine power, the WTGS active power output also linearly 
increases. However, the WTGS active power is slightly smaller than turbine power due to the power loss in the 
induction generator. Figure 8(b) shows that with the increase in turbine power, the rise in WTGS reactive power 
demand is not linear but more rapid or almost exponential. This higher demand for reactive power is because 
more reactive power is needed to magnetize the induction generator as the WTGS active power increases. 


o 
œ 


WPP Active Power Output 
o 
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2 2.5 


Turbine Power 
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Figure 8. Variation of WPP outputs; (a) active power and (b) reactive power 
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Since the WTGS can fulfill some load demands, the power supply (active power) from the distribution 
substation can significantly be reduced as the WTGS active power output increases as shown in Figure 9(a). On 
the other hand, the substation reactive power will increase as the turbine mechanical power (i.e., the WTGS 
active power output) is raised as shown in Figure 9(b). This increase in reactive power due to more reactive 
power is needed for the induction generator magnetization as the WTGS active power increases. 
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Figure 9. Variation of substation outputs (a) active power and (b) reactive power 


5. CONCLUSION 

A method to integrate WTGS into a three-phase unbalanced DSLF analysis has been proposed in this 
paper. The present work extends the single-phase multiple-node model to a three-phase multiple-node model. 
The model is then incorporated into the load flow analysis of a three-phase power distribution network. 
Application of the multiple-node model will only modify the load flow analysis by adding two lines and two 
load buses to the distribution system network where the WTGS is installed. Thus, a standard load flow program 
can be employed to compute the unknown quantities in the DSLF problem formulation. Modification to the 
source code program does not need to be carried out. What needs to be done is to change only the program data 
input due to the additional lines and buses. This paper has used two representative three-phase distribution 
networks (i.e., 19-bus and 25-bus networks) to confirm and validate the proposed method. 
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